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Abstract 



In this paper, we study the sensitivity of CUORE, a bolometric double-beta decay experiment under construction at the 
Laboratori Nazionali del Gran Sasso in Italy. Two approaches to the computation of experimental sensitivity are discussed 
and compared, and the formulas and parameters used in the sensitivity estimates are provided. Assuming a background 
rate of 10 -2 cts/ (keVkgy), we find that, after 5 years of live time, CUORE will have a lcr sensitivity to the neutrinoless 

double-beta decay half-life of T®J 2 (la) = 1.6 X 10 26 y and thus a potential to probe the effective Majorana neutrino mass 

down to 41-95 meV; the sensitivity at 1.64a, which corresponds to 90% C.L., will be 3^,(1.64(7) = 9.5 x 10 25 y. This 
range is compared with the claim of observation of neutrinoless double-beta decay in 76 Ge and the preferred range in 
the neutrino mass parameter space from oscillation results. 

Keywords: neutrino experiment, double-beta decay, sensitivity, bolometer, Poisson statistics 



1. Introduction 

Neutrinoless double-beta decay (Ov/3/3) (see Refs. [2,12, 
Q for recent reviews) is a rare nuclear process hypothesized 
to occur if neutrinos are Majorana particles. In fact, the 
search for 0vf3f3 is currently the only experimentally feasi- 
ble method to establish the Majorana nature of the neu- 
trino. The observation of 0^/3/3 may also probe the abso- 
lute mass of the neutrino and the neutrino mass hierarchy. 
Many experiments, focusing on several different candidate 
decay nuclides and utilizing various detector techniques, 
have sought evidence of this decay 0, 0, ID, 0, H|; next- 
generation detectors are currently under development and 
construction and will begin data taking over the next few 
years. Evidence of 0vf3j3 in 76 Ge has been reported [2, 
ll| but has yet to be confirmed 0, [13, Q El. 

The Cryogenic Underground Observatory for Rare Events 
(CUORE) [13, [13 is designed to search for 0vj3/3 in 130 Tc. 
Crystals made of natural Te02, with an isotopic abun- 
dance of 34.167% of 130 Te [H|, will be operated as bolome- 
ters, serving as source and detector at the same time. Such 
detectors combine excellent energy resolution with low in- 
trinsic background, and they have been operated in stable 
conditions underground for several years 0,0,[Ili. Indi- 
vidual detectors can be produced with masses up to ~ 1 kg, 
allowing for the construction of close-packed large- mass ar- 
rays. Bolometric detectors enable precision measurement 
of the energy spectrum of events inside the crystals, allow- 
ing the search for an excess of events above background 
in a narrow window around the transition energy of the 
isotope of interest. Such a peak constitutes the signature 
of 0v/3/3, and if it is observed, the 0i/[3f3 half-life can be 
determined from the number of observed events. 

The current best limit on 0vB(3 in 130 Te comes from the 
Cuoricino experiment [i, 22, 23 1, which operated 58 crys- 
tals of natural Te02 and 4 enriched Te02 crystals (con- 
taining approximately 11 kg of 130 Te in total) in the Lab- 
oratori Nazionali del Gran Sasso, Italy, from 2003-2008. 
With a total exposure of 19.75 kgy, Cuoricino set a limit 
ofT°fc > 2.8 x 10 24 y (90% C.L. 
of 130 Te. 



jj on the OvfiP half-life 



CUORE, the follow-up experiment to Cuoricino, is cur- 
rently under construction and will exploit the experience 
and results gained from its predecessor. With its 988 de- 
tectors and a mass of - 206 kg of 130 Te, CUORE will be 
larger by more than an order of magnitude. Background 
rates are also expected to be reduced by approximately an 
order of magnitude with respect to Cuoricino. 

In this study, we discuss the sensitivity of CUORE and 
of CUORE-0, the initial phase of CUORE. We start by 
providing the detailed assumptions and formulas for the 
sensitivity estimations. We then review the experimental 
setup and parameters from which the sensitivity values are 
calculated. Finally, we compare the sensitivities with the 
claim of observation of Oz^/3/3 in 76 Ge and the preferred 
range of neutrino masses from oscillation results. 

2. Sensitivity of Double-Beta Decay Experiments 

After introducing some basic 0vf3f3 formulas in Sec. 12. 11 
we will present two possible approaches that can be taken 
to derive an 'experimental sensitivity' that expresses the 
capabilities of an experiment. 

In Sec. 12.21 we will develop a simplified, formula-based 
calculation that uses several basic experimental parame- 
ters (e.g., resolution, background rate, mass) to express 
sensitivity in terms of expected background fluctuations. 
We will refer to this calculation as the "background- 
fluctuation" sensitivity. The background-fluctuation sen- 
sitivity cannot be extended to the ideal zero-background 
case, so we develop an analytical expression for zero- 
background sensitivity in Sec. 12.31 

In Sec. 12.41 we will discuss a procedure to express sensi- 
tivity as the average limit that a particular experiment can 
expect to set in the case that the true 0^/3/3 rate is zero, 
by applying the experiment's analysis tools to a suite of 
Monte Carlo trials. We will refer to this calculation as the 
"average-limit" sensitivity. 

As we will show, the results of the two approaches are 
compatible, although the philosophies of their construction 
differ. 
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2.1. Basic Double-Beta Decay Physics 

Double-beta decay is a second-order weak process, so 
half-lives are typically long: two-neutrino double-beta de- 
cay half-lives are at least of order 10 18 years, while current 
limits on Qvf3/3 half-lives are on the order of 10 24 years or 
greater. With such long half-lives, the radioactive decay 
law can be approximated as 



N(t) ~ N 1 - ln(2) 



t 



T, 



1/2 



(1) 



where Tj/j is the half-life, No is the initial number of 
atoms and N(t) is the number of atoms left after time 
t has passed. 

Assuming that the exchange of a light Majorana neu- 
trino is the dominant 0^/3/3 mechanism, the effective Ma- 
jorana mass of the electron neutrino can be inferred from 
the Oz//3/3 half-life as follows @: 



1/3/3 



1/2 



(2) 



where m e is the electron mass, Fn is a nuclear structure 
factor of merit that includes the nuclear matrix elements 
(NME) and the phase space of the Oj//3/3 transition, and 
T°J 2 is the Ov00 half-life. 

NMEs arc difficult to calculate, and a large range of val- 
ues can be found in the literature, arising from variations 
in the details of the models and the assumptions made. For 
the purpose of this work, we will consider the most recent 
calculations from three different methods: Quasiparticle 
Random Phase Approximation (QRPA) (carried out by 
two different groups: Faessler et at, henceforth denoted 
by QRPA-F, and Suhonen et ai, henceforth denoted by 
QRPA-S), Interacting Shell Model (ISM) and Interacting 
Boson Model (IBM). Fn values and references are shown 
in Tab. [1] These values are calculated using the NMEs 
reported by each group and their suggested phase space 
calculation, taking care to match the values of parameters 
within each model as pointed out in Ref. [24j| . For the 
QRPA calculations, both groups report several values of 
the NMEs depending on the choice of the input parame- 
ters in the model. Therefore, we quote a range of possible 
nuclear factors of merit, taking the maximum and mini- 
mum value reported by each of the two groups. A range is 
also shown for the ISM model to account for the choice of 
different models of short-range correlations. No statistical 
meaning is implied in the use of these ranges. 

2.2. Background- Fluctuation Sensitivity 

The mean value So of the Oz^/3/3 signal, i.e., the expected 
number of Oz^/3/3 decays observed during the live time t is: 



■So 



M ■ Na • a-Tf 
W 



ln(2) 



t 



J l/2 



(3) 



where M is the total active mass, r\ is the stoichiometric 
coefficient of the Oz^/3/3 candidate (i.e., the number of nuclei 



of the candidate Oz//3/3 element per molecule of the active 
mass), W is the molecular weight of the active mass, Na 
is the Avogadro constant, a is the isotopic abundance of 
the candidate Oi//3/3 nuclide and e is the physical detector 
efficiency. 

In Eq. Q, Ty 2 refers to the (unknown) true value of 
the 0^/3/3 half-life, and So is therefore also unknown. In 
the derivation of the background-fluctuation sensitivity, we 
will first determine our sensitivity in terms of a number of 
counts (analogous to So), and then use the form of Eq. §5§ 
to convert to a half-life sensitivity (analogous to Ty 2 ). In 
order to prevent confusion between sensitivities and true 
values, hatted quantities (e.g., T°Y 2 , So) will be used to 
represent the sensitivities corresponding to the unhatted 
true values. 

An experiment can expect to see a background con- 
tribution to the counts acquired in the energy window of 
interest for the Oz^/3/3 signal. In the case of a bolometric ex- 
periment, or indeed any experiment in which the source is 
embedded in the detector (common though not universal 
for Oz//3/3 experiments), we can express the mean number 
of background counts B{SE) in an energy window 5E as 



B(SE) = b-M-SE-t, 



(4) 



where b is the background rate per unit detector mass per 
energy interval (units: cts/(keVkgy)). 

In bolometric experiments, b is independently mea- 
sured, usually by a fit over an energy range much larger 
than the energy window of interest SE. However, the back- 
ground in SE still follows a Poisson distribution with a 
mean value of B(SE). 

An important assumption is implied by the form of 
Eq. namely, that the number of background events 
scales linearly with the absorber mass of the detector. We 
will use this simple model for our background-fluctuation 
sensitivity calculations. However, other cases, most no- 
tably surface contaminations, are in fact possible wherein 
the background might not scale with M. A fully correct 
treatment of an experiment's background would require 
a detailed understanding of the physical distribution of 
the contaminations that are the source of the background, 
used as input for Monte Carlo simulations of the specific 
detector geometry under consideration. 

Throughout the background-fluctuation sensitivity deriva- 
tion, we will use Eq. (O and Eq. Q as analytic expressions 
for the expected numbers of signal and background counts 
assuming a source-equals-detector experimental configura- 
tion, but an analogous estimation is possible for any de- 
tector configuration. 

With the background B(8E) as defined in Eq. (|4]), we 
can calculate the number of counts that would represent 
a positive background fluctuation of a chosen significance 
level. For simplicity, we construct a single-bin counting 
experiment wherein the width of the bin is equal to the 
energy window SE; this way we need to consider only a 
single measured value, sampling a count distribution with 
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Table 1: OufS/3 nuclear factors of merit Fn, as defined in Eq. {2j|, for the candidate Ov/3/3 nuclides discussed in this paper, according to different 
evaluation methods and authors. QRPA: Quasiparticlc Random Phase Approximation; ISM: Interacting Shell Model; IBM: Interacting Boson 
Model. See Se c. I2.1l for details. The phase space values used in calculating Fn values are taken from Table V of Ref. |25| for QRPA-F and 
table 6 of Ref. l26l for the other models. 





Isotope 


QRPA-F H3 


Ovpfi nuclear factor of merit Fn 
(10- 13 y" 1 ) 
QRPA-S [H ISM [H 


IBM-2 [H 


130 Te 

76 Ge 


2.86 - 10.0 
0.88 - 2.60 


2.76 - 7.38 1.86 - 2.91 
0.66 - 1.81 0.33 - 0.50 


6.82 
1.88 



mean B(8E), and we can decouple the sensitivity calcu- 
lation from the specific analysis approach used by the ex- 
periment. 

We can now define our background-fluctuation sensitiv- 
ity: it is the smallest mean signal So that is greater than or 
equal to a background fluctuation of a chosen significance 
level. It is common in this kind of study [3l| to express the 
significance level in terms of Gaussian standard deviations 
from the background value. We will speak in general terms 
of 'sensitivity at n a \ where n a is the desired significance 
level in terms of number of Gaussian a. Thus, if B(8E) is 
large enough that the background count distribution can 
be considered to be Gaussian, the desired value of S(SE) is 
determined by setting the following requirement in terms 
of cr = y/B(5E): 



S(SE) = So ■ f(SE) = n a ■ y/B(8E), 



(5) 



where f(SE) is the fraction of signal events that fall in 
the energy window cut 5E around the Q-value. The inclu- 
sion of f{SE) arises from our construction of a single-bin 
counting experiment; it serves as a simple estimate of the 
analysis efficiency. 

For a Gaussian signal (i.e., Gaussian-distributed in en- 
ergy around the Q-value), the signal fraction f(5E) is 



f(SE) = erf 



/ SE 
\~AE 



(0) 



where AE is the detector FWHM energy resolution. The 
value of SE can be chosen to maximize the S(5E)-to- 
y / 'B(SE) ratio in the energy window of interest, which 
in turn optimizes the sensitivity criterion expressed by 
Eq. ([5]); this optimal choice corresponds to SE w 1.2AE. 
It is, however, common to take SE = AE. In this case, the 
sensitivity differs by less than 1% from the one calculated 
at the optimal cut. 

By using the expressions for So and B(SE) from Eq. ^ 
and ((U), we obtain the Gaussian-regime expression for the 
background-fluctuation sensitivity of experiments in 
the following form: 



ln(2) N A ■ a ■ 77 ■ e M ■ t 
~rL W V b-SE 



f(SE). (7) 



This equation is extremely useful in evaluating the ex- 
pected performances of prospective experiments, as it an- 
alytically links the experimental sensitivity with the de- 
tector parameters. Aside from the inclusion of the signal 
fraction, it is similar to the familiar 'factor of merit' ex- 
pression used within the Oz//3/3 experimental community. 

For small numbers of observed events, the Gaussian 
approximation of Eq. and Eq. (O does not provide the 
correct probability coverage, and therefore the meaning of 
the significance level is not preserved. In fact, the Gaus- 
sian approximation for the distribution of the number of 
observed counts becomes invalid when the expected num- 
ber of background counts is small; if B{SE) is less than 
~ 24 counts, the Gaussian calculation of a ler sensitiv- 
ity will differ from its Poissonian counterpart (developed 
below) by 10% or more. 

Although the Gaussian limit will possibly still be suf- 
ficient for CUORE (see SecSJ), a more careful calculation 
might be necessary in the case of a lower background or 
smaller exposure, or for more sensitive experiments in the 
future. We therefore compute the sensitivity by assuming 
a Poisson distribution of the background counts. 

In terms of Poisson-distributed variables, the concept 
expressed by Eq. §5§ becomes [32 



E 



PB (k) = a, (8) 

k=S(SE)+B{SE) 

where a is the Poisson integrated probability that the 
background distribution alone will cause a given exper- 
iment to observe a total number of counts larger than 
S{8E) + B(SE). As written, Eq. (JSJ can be solved only for 
certain values of a because the left-hand side is a discrete 
sum. To obtain a continuous equation that preserves the 
Poisson interpretation of Eq. © , we exploit the fact that 
the (discrete) left-hand side of Eq. @ coincides with the 
(continuous) normalized lower incomplete gamma function 
P(a, x) (see page 260 of Ref. [33| for details): 

P(S(SE)+B(SE) 1 B(SE))=a. (9) 

The computation of So from Eq. @, for given values of 
B(SE) and a, can be achieved numerically. Once So is 
computed in this way, the corresponding Poisson-rcgimc 
background-fluctuation sensitivity to the half-life T®Y 2 for 
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ncutrinolcss doublc-bcta decay is simply calculated by re- 
versing Eq. ([3]). 

For the remainder of this paper, we will use the Poisson- 
regime calculation based on Eq. © to evaluate our 
background-fluctuation sensitivity. So that we can con- 
tinue to indicate the significance level with the familiar 
n a notation instead of the less-intuitive a, however, we 
will label our sensitivities with the n a corresponding to a 
Gaussian upper-tail probability of a (for example, we will 
call a background-fluctuation sensitivity calculated with 
a = 0.159 in Eq. © a 'la sensitivity'). 



2.3. Analytical Expression for Zero- Background Sensitiv- 
ity 

It is meaningless to define sensitivity in terms of back- 
ground fluctuations when B(SE) = 0; therefore, the background 
fluctuation sensitivity calculation cannot be extended to 
the ideal 'zero-background' case. If we wish to develop a 
formula-based, analysis-decoupled zero-background sensi- 
tivity calculation, we can still construct a single-bin count- 
ing experiment in the same way as we did for the background- 
fluctuation sensitivity; however, we must adopt a new method 
of constructing our sensitivity parameter. 

To construct the zero-background sensitivity, we choose 
to follow the Bayesian limit-setting procedure. Instead 
of comparing the mean signal value S(5E) to the mean 
background value B(5E), we are now obliged to consider 
S m ax(SE), the upper limit on S(5E) in the case that the 
experiment observes zero counts (i.e., no background or 
signal) in SE during its live time. S max (6E) can be eval- 
uated using a Bayesian calculation with a flat signal prior 
(see Eq. (32.32)-(32.34) of Ref. (H): 



2.4- Average-Limit Sensitivity 

The average-limit sensitivity calculation is a Monte- 
Carlo-based procedure constructed in a similar manner as 
the zero-background sensitivity presented in the previous 
section. Following what we have done in the method 
works as follows: 



1. Generate a large number of toy Monte Carlo spectra 
assuming zero Qvj3f3 signal in the fit window (much 
wider than the 5E = AE window used for the 
background-fluctuation sensitivity, in order to utilize 
the available shape information in the fit). 

2. For each Monte Carlo spectrum, perform a binned 
maximum likelihood fit to the spectrum and extract 
the associated Bayesian limit with a flat signal prior 
by integrating the posterior probability density (the 
same analysis technique used in 0, [23). 

3. Construct the distribution of the limits calculated 
from the Monte Carlo spectrum, and determine its 
median. 



rS ma .. 

Js=o 



(SE) 



Ps(0)dS 



rS ma: 

Js=o 



(SE) 



S°e~ s dS C.L. 



J™ =o Ps{0)dS 



f™ S°eSdS 



100 



(10) 



where ps{k) is the Poisson distribution p M (fc) with mean 
(i = S and the credibility level C.L. is expressed as a per- 
cent. Eq. (fTU|) can be solved analytically for S max (5E): 

S max (SE) = S max ■ f(5E) = - ln(l - ^), (11) 

where S max is the inferred upper limit on So- Using S max 
in place of Sq in Eq. ©, we obtain 



T°; 2 (C, 



L.) = - 



ln(2) N A -a-r) 



Hi « 



W 



-M-t-f(SE). (12) 



For practical purposes, this background-free approx- 
imation becomes valid when the expected value of the 
background is of the order of unity, B(8E) < 1 count. 
It should be stressed that, because the zero-background 
sensitivity is (by necessity) constructed differently than 
the background-fluctuation sensitivity, the interpretations 
of the two do not entirely coincide. 



The average-limit sensitivity method is, in a way, more 
powerful than the analytical background-fluctuation method 
because it can in principle take into account subtle and 
detector-dependent experimental effects, which can be dif- 
ficult or sometimes impossible to model with analytical 
formulas. However, because the average-limit approach re- 
lies on analysis of statistical ensembles, it lacks the great 
advantages of clarity and simplicity offered by straight- 
forward formulas. It is clear that the two methods must 
be (and indeed are, as shown later) essentially equivalent 
given the same input parameters, though a minor system- 
atic difference arises because the probability distribution 
of the limits is not symmetric and the median found with 
the MC does not coincide with the S(SE) computed with 
Eq. @. 

For a completed experiment like Cuoricino, the experi- 
mental parameters (e.g., background rate(s) and shape(s), 
resolution(s), exposure) used as inputs to the Monte Carlo 
in step [1] are the real parameters that have been directly 
measured by the experiment. The average- limit sensitiv- 
ity is meaningful for a completed experiment that has not 
seen evidence of a signal because it provides an under- 
standing of how 'lucky' the experiment was in the limit 
it was able to set. To adapt the approach for an upcom- 
ing experiment, it is of course necessary to instead use the 
expected experimental parameters to generate the Monte 
Carlo spectra in step[TJ Calculating the average-limit sen- 
sitivity in this way allows for the direct comparison of an 
upcoming experiment with previously reported experimen- 
tal limits. The average-limit sensitivity is often also the 
value 0^/3/3 physicists have in mind when they consider 
the meaning of sensitivity; for example, the GERDA ex- 
periment reports a sensitivity calculated in essentially this 
although they choose to report the mean ex- 



manner 



24, 



pected limit instead of the median. 
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3. Validation of the Methods with Cuoricino 

Cuoricino [35| achieved the greatest sensitivity of any 
bolometric 0vf3f3 experiment to date and served as a pro- 
totype for the CUORE experiment. Cuoricino took data 
from 2003 to 2008 in the underground facilities of the Lab- 
oratori Nazionali del Gran Sasso (LNGS), Italy. 

The Cuoricino detector consisted of 62 Te02 bolome- 
ters with a total mass of 40.7 kg. The majority of the 
detectors had a size of 5 x 5 x 5 cm 3 (790 g) and consisted 
of natural Te02- The average FWHM resolution for these 
crystals was 6.3 ± 2.5 keV at 2615 keV Q, the nearest 
strong peak to the 0vf3f3 transition energy. Their physi- 
cal efficiency, which is mostly due to the geometrical effect 
of beta particles escaping the detector and radiative pro- 
cesses, has been estimated to be Ephys = 0.874 ±0.011 Li- 
The full details of the crystal types present in the detector 
array can be found in |4|. 

The most recent Cuoricino limit was published along- 
side an average-limit sensitivity. This sensitivity was eval- 
uated as the median of the distribution of 90% C.L. limits 
extracted from toy Monte Carlo simulations that used the 
measured detector parameters as inputs, and it was deter- 
mined to be 2^ 2 (90% C.L.) = 2.6 x 10 24 y. 

Because of the different crystal types present in Cuori- 
cino, if we wish to calculate a background-fluctuation sen- 
sitivity for Cuoricino to compare with this average-limit 
sensitivity, we need to slightly adjust the background- 
fluctuation calculation presented in Sec. 12.21 to accommo- 
date different parameter values for the different crystal 
types. Cuoricino can be considered as the sum of virtual 
detectors, each representing one of the crystal types during 
one of two major data-taking periods, called Runs. The 
detectors' total exposures, background rates after event 
selection, physical efficiencies, and average resolutions are 
reported in Ref. Q , subdivided by crystal type and Run as 
appropriate. Therefore we can use these reported values 
to calculate both our expected signal S(SE) and expected 
background B(SE) as sums of the contributions from these 
virtual detectors, then follow the Poisson-regime background- 
fluctuation sensitivity procedure. If we wish our background- 
fluctuation sensitivity to be quantitatively comparable to 
a 90% C.L. average-limit sensitivity, we must choose to 
calculate the background-fluctuation sensitivity at 1.64a 

{a = 0.051); indeed, doing so yields 2^(1.64a) = 2.6 x 
10 24 y, in perfect agreement with the average-limit sensi- 
tivity. 

Following previously established convention for past 
bolometric experiments (l7l . 136^ . we choose to report 
background-fluctuation sensitivities at la (a = 0.159) for 
upcoming experiments. For the purpose of illustration, 
the corresponding background-fluctuation sensitivity for 
Cuoricino would be T°f 2 (la) = 4.2 x 10 24 y. 

Although upcoming CUORE-family experiments have 
historically shown la background-fluctuation sensitivities, 
which quantitatively roughly coincide with 68% C.L. average- 



limit sensitivities, other upcoming Oz//3/3 experiments com- 
monly report 90% C.L. sensitivities. To prevent confusion 
between our sensitivity approach and that commonly used 
by other Ov/3/3 experiments, it is instructive to compare 
1.64(7 background-fluctuation sensitivities to 90% C.L. 
average-limit sensitivities for both CUORE and CUORE- 
0; this comparison appears in Sec. [4j 

4. CUORE sensitivity 

CUORE will consist of an array of 988 Te02 cubic 
detectors, similar to the 5x5x5 cm 3 Cuoricino crys- 
tals described above. The total mass of the detectors 
will be 741 kg. The detectors will be arranged in 19 in- 
dividual towers and operated at ~ 10 mK in the Gran 
Sasso underground laboratory. The expected energy res- 
olution FWHM of the CUORE detectors is AE « 5 keV 
at the QvBB transition energy, or Q-value (~ 2528 keV 
fo: 



130r 



Te 37, 38 



391]). This resolution represents an im- 
provement over that seen in Cuoricino and has already 
been achieved in tests performed in the CUORE R&D fa- 
cility at LNGS. CUORE is expected to accumulate data 
for about 5 years of total live time. The experiment is cur- 
rently being constructed and first data-taking is scheduled 
for 2014. 

The CUORE collaboration plans to operate a single 
CUORE-like tower in the former Cuoricino cryostat, start- 
ing in late 2011. This configuration, named CUORE-0, 
will validate the assembly procedure and the readiness of 
the background reduction measures. The experimental pa- 
rameters of CUORE-0 and CUORE that are used in the 
sensitivity calculations are summarized in Tab. [2] 

The background rate is the most critical parameter to 
assess before the calculation of the sensitivity can be car- 
ried out. 

In Cuoricino, the average background counting rate 
in the region of interest (ROI) for Oi//3/3 decay, namely, 
a region centered at the Q-value and 60 keV wide, was 
0.161 ±0.006 cts/(keVkgy) for the 5 x 5 x 5 cm 3 crystalfl 
An analysis of the background sources responsible for the 
flat background in the ROI has been performed on a par- 
tial set of statistics [It], , following the technique and 
the model developed for the MiDBD experiment |40| . The 
result of this analysis was the identification of three main 
contributions: 30 ± 10% of the measured fiat background 
in the ROI is due to multi-Compton events due to the 
2615 keV gamma ray from the decay chain of 232 Th from 
the contamination of the cryostat shields; 10 ±5% is due to 
surface contamination of the Te02 crystals with 238 U and 
232 Th (primarily degraded alphas from these chains); and 
50 ± 20% is ascribed to similar surface contamination of in- 
ert materials surrounding the crystals, most likely copper 



6 This is the background rate measured when operating the array 
in anticoincidence; this evaluation is extracted from the Of/3/3 best 
fit p| and corrected for the instrumental efficiency to give the real 
rate. 
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Table 2: Values used in the estimation of the sensitivity of CUORE-0 and CUORE. Symbols arc defined in Eq. (0, Eq. {4j, and Eq. JSJl. See 
Sec.[4]for a discussion of the background values. 



Experiment 


a 

(%) 


V 


e 

(%) 


W 

(g/mol) 


M 

(kg) 


AE 

(kcV) 


/(AS) 

(%) 


b 

(cts/(kcVkgy)) 


CUORE-0 
CUORE 


34.167 
34.167 


1 
1 


87.4 
87.4 


159.6 
159.6 


39 
741 


5 
5 


76 
76 


0.05 
0.01 





(other sources that could contribute are muons [41| and 
neutrons, but simulations indicate that these have only a 
minor effect). 

On the basis of this result, the R&D for CUORE has 
pursued two major complementary avenues: one, the re- 
duction of surface contamination, and two, the creation 
of an experimental setup in which potential background 
contributions arc minimized by the selection of extremely 
radio-pure construction materials and the use of highly 
efficient shields. The latter activity is based mainly on 
standard procedures (material selection with HPGe spec- 
troscopy, underground storage to avoid activation, eval- 
uation of the background suppression efficiencies of the 
shields on the basis of Monte Carlo simulations (42[, etc.). 
However, the required surface contamination levels are ex- 
tremely low, on the order of 1-10 nBq/cm 2 , nearly un- 
detectable with any standard technique used in surface 
analysis. In most cases, only bolometric detectors are suf- 
ficiently sensitive; at this time, our understanding of these 
contaminations comes only from the statistics-limited data 
sets collected by small test detectors constructed from 
CUORE materials (see Ref. 43( for the contract require- 
ments on and measurements of the contamination levels of 
the crystals). 

A detailed analysis of the background mitigation effort 
and its extrapolation to the CUORE and CUORE-0 back- 
ground is out of the scope of the present paper. A full 
account of the performed measurements, analysis, and re- 
sults is being prepared and will be published soon. Here, 
to justify the expected background rates that will be used 
for the sensitivity estimations, we offer a brief summary, 
allowing us to perform a simple scaling to obtain the range 
into which we expect the CUORE-0 background rate to 
fall and support the conclusion that CUORE will meet its 
design background specification. 

CUORE crystals are produced following a controlled 
protocol [44[ that is able to ensure a bulk contamination 
level lower than 3 x 10 -12 g/g in both 238 U and 232 Th. A 
more rigorous surface-treatment technique than that used 
for the Cuoricino crystals was developed; when studied 
with a small array of bolometric detectors, it proved to 
be able to reduce the surface contamination of Cuoricino 
crystals re-treated with this method by approximately a 
factor of 4 [45j]. The technique has now been adopted 
and applied in the production of the CUORE crystals, 
and bolometric tests have already proven its efficacy (43|. 
A preliminary evaluation of the surface contaminations of 



the final CUORE crystals [44| indicated a lower limit on 
the reduction with respect to the contamination seen in 
Cuoricino of a factor of 2; the measurement was statistics- 
limited, so the true reduction factor may be greater. 

In Cuoricino, a large fraction of the 0^/3/3 background 
was identified as due to surface contamination of the cop- 
per — the only significant material surrounding the detec- 
tors, which are mounted in vacuum. Unfortunately, the 
signature of the surface contamination of the copper is ex- 
tremely weak when compared to other contributions, as 
the background ascribed to the copper contamination is 
a flat continuum that can be easily observed only in the 
peakless 3-4 MeV region of the spectrum [4(J |45[ . Exten- 
sive efforts have been dedicated to the study of different 
treatment procedures able to reduce the copper surface 
contamination; in the end, a technique that proved to be 
capable of reducing the copper surface contamination by 
at least a factor of 2 as compared with that observed in 
Cuoricino has been selected by the collaboration as the 
baseline for the CUORE copper treatment. 

Based on the above-reported considerations, we may 
define a conservative case wherein we assume that the spe- 
cific contaminations of the CUORE copper and crystals 
have both been reduced by a factor of 2 relative to Cuori- 
cino. There is a good chance that the reduction factors 
are much higher, but this cannot be confirmed at present 
due to the limited statistics of our measurements. Only 
CUORE-0 will ultimately be able to measure the true level 
of radiopurity achieved with the chosen surface treatment. 

CUORE-0 will consist of CUORE crystals mounted in 
CUORE-style frames as a single tower. Because of this 
geometry, which is similar to that of Cuoricino, the con- 
tamination reduction factors reported above scale almost 
directly to the background we expect to observe in the 
ROI. The total amount of copper facing the crystals will 
be only slightly reduced with respect to Cuoricino, but 
its surface will be treated with the new procedure studied 
for CUORE. CUORE-0 will be assembled in the Cuori- 
cino cryostat, so the gamma background from contamina- 
tion in the cryostat shields will remain approximately the 
same as in Cuoricino. We consider that the irreducible 
background for CUORE-0 comes from the 2615 kcV 208 T1 
line due to 232 Th contaminations in the cryostat, in the 
case that all other background sources (i.e., surface con- 
taminations) have been rendered negligible; this would 
imply a lower limit of ~ 0.05 cts/(keVkgy) on the ex- 
pected background in CUORE-0. Similarly, an upper limit 
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Figure 1: CUORE-0 background- fluctuation sensitivity at lcr for 
two different values of the background rate in the region of inter- 
est, 0.05 cts/(keVkgy) (solid line) and 0.11 cts/(keVkgy) (dotted 
line), representing the range into which the CUORE-0 background 
is expected to fall. 



Figure 2: la expected background-fluctuation sensitivities for the 
CUORE-0 (dotted line) and CUORE (solid line) experiments, cal- 
culated from Eq. (0 and Eq. with the experimental parameters 
shown in Tab. [2] The Cuoricino la sensitivity calculation (dashed 
line) is discussed in Sec. \3\ 



of 0.11 cts/(keVkgy) follows from scaling the Cuoricino 
background in the conservative case, described above, of a 
factor of 2 improvement in crystal and copper contamina- 
tion. 

A plot of the expected lcr background-fluctuation sen- 
sitivity of CUORE-0 as a function of live time in these two 
bounding cases is shown in Fig.[T] Tab. [3] provides a quan- 
titative comparison between la background-fluctuation sen- 
sitivities (as shown in Fig. [1]), 1.64cr background- fluctuation 
sensitivities, and 90% C.L. average-limit sensitivities for 
CUORE-0 at several representative live times. The antic- 
ipated total live time of CUORE-0 is approximately two 
years; for this live time at the 0.05 cts/(keVkgy) back- 
ground level, B(SE) ~ 20 cts, meaning that the Poisson- 
rcgime calculation is really necessary in this case because it 
differs from the Gaussian-regime approximation by > 10% 
(see Sec. [22]). 

CUORE, in addition to the new crystals and frames 
already present in CUORE-0, will be assembled as a 19- 
tower array in a newly constructed cryostat. The change 
in detector geometry will have two effects. First, the large, 
close-packed array will enable significant improvement in 
the anticoincidence analysis, further reducing crystal-related 
backgrounds. Second, the fraction of the total crystal sur- 
face area facing the outer copper shields will be reduced by 
approximately a factor of 3. In addition to these consider- 
ations, the new cryostat will contain thicker lead shielding 
and be constructed of cleaner material, which should result 
in a gamma background approximately an order of magni- 
tude lower than that in the Cuoricino cryostat. Based 
on the above considerations and the Cuoricino results, 
CUORE is expected to achieve its design background value 
of 0.01 cts/(keVkgy). A comprehensive Monte Carlo sim- 
ulation that includes the most recent background measure- 
ments is currently ongoing. 

An overview of the lcr background-fluctuation sensitiv- 



ities of the Cuoricino, CUORE-0, and CUORE Te0 2 bolo- 
metric experiments is shown in Fig. [2] The Cuoricino lcr 
sensitivity calculated in Sec. [3] is shown for reference. A 
lcr half-life sensitivity close to 10 25 years is expected from 
2 years' live time of CUORE-0. Once CUORE starts data- 
taking, another order of magnitude improvement in sensi- 
tivity is expected in another two years. 

A plot of the CUORE experiment's sensitivity as a 
function of the live time and exposure is shown in Fig. [3j 
Tab. 2] provides a quantitative comparison between lcr 
background-fluctuation sensitivities (as shown in Fig. [3]), 
1.64cr background-fluctuation sensitivities, and 90% C.L. 
average-limit sensitivities for CUORE at several repre- 
sentative live times. The anticipated total live time of 
CUORE is approximately five years; for this live time at 
the design goal background level, B(5E) ~ 190 cts, mean- 
ing that the Gaussian approximation would still be valid 
in this case. The sensitivity values we show in this pa- 
per nevertheless differ from those previously reported by 
the experiment 1(1 17 1, but this ~ 25% difference can be 



attributed to the inclusion of the signal fraction f(6E), 
which has not previously been considered. 

As mentioned previously, estimates of the CUORE back- 
ground arc currently based on measured limits, not mea- 
sured values. While there are promising indications that it 
may perform even better than its design value of 
0.01 cts/ (keV kg y) , it is also likely that background rates 
of 0.001 cts/(keVkgy) or below cannot be reached with 
the present technology. Even so, R&D activities are al- 
ready underway pursuing ideas for further reduction of 
the background in a possible future experiment. Tech- 
niques for active background rejection are being investi- 
gated 0, 47 1) that could provide substantial reduction 
of the background. Sensitivities for a scenario with 0.001 
cts/(keVkgy) in a CUORE-likc experiment are given in 
Fig. [3] and Tab.H 



8 



Table 3: Background-fluctuation half-life sensitivities at lcr for CUORE-0 under different background estimations after one, two, and four 
years of live time. The boldcd column corresponds to the approximate anticipated total live time of two years. 1.64<r background- fluctuation 
sensitivities and 90% C.L. average-limit sensitivities, in italics, are also provided to illustrate the similarity of the two values. 
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(kcV) 


(sig./conf. level) 
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4y 


0.11 


5 


la 


0.45 


0.66 


0.95 






1.64a 


0.28 


0.40 


0.58 






90% C.L. 


0.29 


0.41 


0.59 


0.05 


5 


la 


0.64 


0.94 


1.4 






1.64a 


0.39 


0.58 


0.84 






90% C.L. 


0.39 


0.59 


0.83 



Table 4: Background-fluctuation half-life sensitivities at lcr for CUORE after two, five, and ten years of live time. The bolded column 
corresponds to the approximate anticipated total live time of five years. The sensitivities are reported for the design goal background level, 
as well as for an order-of-magnitude improvement over the design goal. 1.64<r background-fluctuation sensitivities and 90% C.L. average-limit 
sensitivities, in italics, are also provided to illustrate the similarity of the two values. 
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5. Comparison with the claim in 76 Ge 
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Figure 3: Baseline expected background-fluctuation sensitivity of 
the CUORE experiment at lcr (solid line). The sensitivity for an 
order-of-magnitude improvement over the baseline background is also 
shown (dotted line). 



It is interesting to compare the CUORE-0 and CUORE 
sensitivities with the claim for observation of OV/3/3 in 76 Ge 0, 
Tot 11 1 . The authors of this claim have reported several 



different values for the half-life of 76 Ge, depending upon 
the specifics of the analysis; the longest of these, and 
thus the one requiring the greatest sensitivity to probe, 



is T°; 2 ( 76 Ge) 
follows that 



2.23 



+0.44 

-0.31 



10 25 y From Eq. ©, it 



-i0v /130 



1/2 



F N ( 7b Ge) 
F N ( 130 Te) 



(76 



1/2 



( 7fa Ge). 



However, because of the wide spread in Fn calculations 
(see Tab. [1]), directly using this equation to estimate the 
expected half-life for the 0i/(3[3 of 130 Te can be mislead- 
ing. A method of treating NME uncertainties based on 
the QRPA-F calculations is suggested, and shown to be 
roughly consistent with the QRPA-S and ISM calculations, 
in [25J. Following this method, the expected lcr range of 
T°J 2 ( 130 Te) is (5.1 - 7.1) x 10 24 y (for the best-fit value of 
the 76 Ge claim) or (4.2 - 8.9) x 10 24 y (when the lcr uncer- 
tainty on the 76 Ge claim is also included, slightly shifting 
the central value of the log 10 [T^ ( 76 Ge)] range as done 
in pBl ] so that the errors become symmetric). 

The mathematical framework of the background- 
fluctuation sensitivity calculation can be inverted to de- 
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Figure 4: Significance level at which CUORE-0 can observe a signal 
corresponding to the 76 Ge claim, assuming the best expected back- 
ground of 0.05 cts/ (kev kg y ) . The inner band assumes the best- fit 
value of the 76 Ge claim, and its width arises from the la range of 
QRPA-F NMEs calculated in [25|. The outer band accounts for the 
ltr uncertainty on the 76 Ge claim in addition to the range of NMEs. 



termine the magnitude of the mean signal in terms of n a 
that an assumed 'true' half-life value will produce in an ex- 
periment. Fig. 2] shows the n a significance level at which 
CUORE-0 can probe the 76 Ge claim as it accrues statistics 
over its anticipated live time. The width of the inner band 
corresponds to the la range of NMEs determined in (ijj . 
while the outer band includes the uncertainty on the claim 
as well; i.e., each band is bounded by curves corresponding 
to the maximum and minimum Ty 2 ( 130 Te) of its respec- 
tive range, as given above. As can be deduced from the 
plot, CUORE-0 will achieve at least a la sensitivity to 
any signal within the expected la range of T®J 2 ( 130 Te) 
within two years. By combining data from CUORE-0 and 
Cuoricino, the claim could be verified in a shorter time 
with higher sensitivity. 

Thanks to the increased size and lower background, if 
the 130 Te Of/3/3 half-life indeed falls in the la range im- 
plied by the claim in 76 Ge, CUORE will already be able 
to achieve a 5<r expected signal above background within 
about six months. 



6. Conclusions 

In recent years, experimenters have made great strides 
in the search for neutrinoless double-beta decay, a discov- 
ery which would establish the Majorana nature of the neu- 
trino and have far-reaching ramifications in physics. Next- 
generation Of/3/3 experiments like CUORE have two pri- 
mary goals: to test the claim of observation of Of/3/3 in 
76 Ge, and to begin to probe effective neutrino masses of 
m /3/3 < 50 meV (commonly referred to as the 'inverted 
hierarchy region' of the neutrino mass phase space). We 
have investigated the expected performance of CUORE, 
allowing evaluation of its ability to meet these two goals. 



In Sec.0 we developed two different approaches to cal- 
culating experimental sensitivity: the background-fluctuation 
sensitivity and the average-limit sensitivity. The background- 
fluctuation sensitivity parametrizes the signal that the ex- 
periment is capable of observing in terms of the expected 
background fluctuations, while the average- limit sensitiv- 
ity is the average limit that the experiment expects to 
set in the case that there is no signal to find. Although 
the average-limit sensitivity is more directly comparable 
to previously reported limits by construction, we prefer to 
evaluate upcoming experiments in terms of the background- 
fluctuation sensitivity because the goal of Of/3/3 experi- 
ments is to discover and measure neutrinoless double beta 
decay, not merely set a limit. In fact, the two methods pro- 
duce quantitatively similar results, so it is not misleading 
to consider the background-fluctuation sensitivity as an 
approximation of the average- limit sensitivity if the signif- 
icance/credibility levels of the two methods are properly 
chosen to coincide. 

Tab.[5]contains a summary of la background-fluctuation 
sensitivities to the neutrino Majorana mass according to 
different NME calculations, assuming that the exchange 
of a light Majorana neutrino is the dominant Of /3/3 mech- 
anism, as discussed in Sec. 12.11 These values are consid- 
ered the official sensitivity values for CUORE-family ex- 
periments. During its run, CUORE will fully explore the 
130 Te Of/3/3 half-life range corresponding to the the claim 
of observation of Of/3/3 in 76 Ge. 

For illustrative purposes, Tab. [5] also shows the limiting 
"zero-background" case for both CUORE-0 and CUORE. 
The calculation is performed at 68% C.L. so that the val- 
ues can be considered as zero-background extrapolations 
of the finite-background la background-fluctuation sen- 
sitivities. CUORE-0 and CUORE will both have suffi- 
ciently good resolution that the signal fraction may be 
omitted from Eq. (|T2"j) for the calculation. As discussed 
in Sec. 12.31 the zero-background approximation applies 
when B(SE) < 1 count; we can determine the background 
rate that each experiment would have to achieve to ful- 
fill this requirement from Eq. ([?]), assuming a window of 
5E = 2.5AE (large enough that including f(SE) would 
not change the values reported in Tab. [5]). CUORE-0 
would require b < 1.0 x 10~ 3 cts/(keVkgy); CUORE 
would require b < 2.2 x 10~ 5 cts/(keVkgy), nearly three 
orders of magnitude better than the baseline background 
rate. 

In Fig. [5J the expected sensitivity of CUORE is com- 
pared with the preferred values of the neutrino mass pa- 
rameters obtained from neutrino oscillation experiments. 
The sensitivity of CUORE will allow the investigation of 
the upper region of the effective Majorana neutrino mass 
phase space corresponding to the inverted hierarchy of neu- 
trino masses. 



10 



Table 5: Summary table of expected parameters and la background-fluctuation sensitivity in half-life and effective Majorana neutrino mass. 
The different values of map depend on the different NME calculations; see Sec. 12, ll and Tab. ^ Zero-background sensitivities, in italics, are 
also provided as an estimation of the ideal limit of the detectors' capabilities; they are presented at 68% C.L. so that they can be considered 
as approximate extrapolations of the la background-fluctuation sensitivities. 
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Figure 5: The Cuoricino result and the expected CUORE la background-fluctuation sensitivity overlaid on plots that show the bands 
preferred by neutrino oscillation data (inner bands represent best-fit data; outer bands represent data allowing 3<r errors) [4H , Both normal 
(Am| 3 > 0) and inverted (Am| 3 < 0) neutrino mass hierarchies are shown, (a) The coordinate plane represents the parameter space of 
map and mughtesty following the plotting convention of [4H , (b) The coordinate plane represents the parameter space of mpp and Sm;, 
following the plotting convention of [49| . The widths of the Cuoricino and CUORE bands are determined by the maximum and minimum 
values of man obtained from the four NME calculations considered in this work. 
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